A systematic pipeline of protein structure selection for computer-aided drug discovery: A case study on T790M/L858R mutant EGFR structures.

Virtual screening (VS) is a routine method to evaluate chemical libraries for lead identification. Therefore, the selection of appropriate protein structures for VS is an essential prerequisite to identify true actives during docking. But the presence of several crystal structures of the same protein makes it difficult to select one or few structures rationally for screening. Therefore, a computational prioritization protocol has been developed for shortlisting crystal structures that identify true active molecules with better efficiency. As identification of small-molecule inhibitors is an important clinical requirement for the T790M/L858R (TMLR) EGFR mutant, it has been selected as a case study. The approach involves cross-docking of 21 co-crystal ligands with all the structures of the same protein to select structures that dock non-native ligands with lower RMSD. The cross docking performance was then correlated with ligand similarity and binding-site conformational similarity. Eventually, structures were shortlisted by integrating cross-docking performance, and ligand and binding-site similarity. Thereafter, binding pose metadynamics was employed to identify structures having stable co-crystal ligands in their respective binding pockets. Finally, different enrichment metrics like BEDROC, RIE, AUAC, and EF1% were evaluated leading to the identification of five TMLR structures (5HCX, 5CAN, 5CAP, 5CAS, and 5CAO). These structures docked a number of non-native ligands with low RMSD, contain structurally dissimilar ligands, have conformationally dissimilar binding sites, harbor stable co-crystal ligands, and also identify true actives early. The present approach can be implemented for shortlisting protein targets of any other important therapeutic kinases.

Benchmarking Refined and Unrefined AlphaFold2 Structures for Hit Discovery.

The recently developed AlphaFold2 (AF2) algorithm predicts proteins' 3D structures from amino acid sequences. The open AlphaFold protein structure database covers the complete human proteome. Using an industry-leading molecular docking method (Glide), we investigated the virtual screening performance of 37 common drug targets, each with an AF2 structure and known holo and apo structures from the DUD-E data set. In a subset of 27 targets where the AF2 structures are suitable for refinement, the AF2 structures show comparable early enrichment of known active compounds (avg. EF 1%: 13.0) to apo structures (avg. EF 1%: 11.4) while falling behind early enrichment of the holo structures (avg. EF 1%: 24.2). With an induced-fit protocol (IFD-MD), we can refine the AF2 structures using an aligned known binding ligand as the template to improve the performance in structure-based virtual screening (avg. EF 1%: 18.9). Glide-generated docking poses of known binding ligands can also be used as templates for IFD-MD, achieving similar improvements (avg. EF 1% 18.0). Thus, with proper preparation and refinement, AF2 structures show considerable promise for in silico hit identification.

Computational identification of natural product inhibitors against EGFR double mutant (T790M/L858R) by integrating ADMET, machine learning, molecular docking and a dynamics approach.

Double mutated epidermal growth factor receptor is a clinically important target for addressing drug resistance in lung cancer treatment. Therefore, discovering new inhibitors against the T790M/L858R (TMLR) resistant mutation is ongoing globally. In the present study, nearly 150 000 molecules from various natural product libraries were screened by employing different ligand and structure-based techniques. Initially, the library was filtered to identify drug-like molecules, which were subjected to a machine learning based classification model to identify molecules with a higher probability of having anti-cancer activity. Simultaneously, rules for constrained docking were derived from three-dimensional protein-ligand complexes and thereafter, constrained docking was undertaken, followed by HYDE binding affinity assessment. As a result, three molecules that resemble interactions similar to the co-crystallized complex were selected and subjected to 100 ns molecular dynamics simulation for stability analysis. The interaction analysis for the 100 ns simulation period showed that the leads exhibit the conserved hydrogen bond interaction with Gln791 and Met793 as in the co-crystal ligand. Also, the study indicated that Y-shaped molecules are preferred in the binding pocket as it enables them to occupy both pockets. The MMGBSA binding energy calculations revealed that the molecules have comparable binding energy to the native ligand. The present study has enabled the identification of a few ADMET adherent leads from natural products that exhibit the potential to inhibit the double mutated drug-resistant EGFR.

An electron transfer competent structural ensemble of membrane-bound cytochrome P450 1A1 and cytochrome P450 oxidoreductase.

Cytochrome P450 (CYP) heme monooxygenases require two electrons for their catalytic cycle. For mammalian microsomal CYPs, key enzymes for xenobiotic metabolism and steroidogenesis and important drug targets and biocatalysts, the electrons are transferred by NADPH-cytochrome P450 oxidoreductase (CPR). No structure of a mammalian CYP-CPR complex has been solved experimentally, hindering understanding of the determinants of electron transfer (ET), which is often rate-limiting for CYP reactions. Here, we investigated the interactions between membrane-bound CYP 1A1, an antitumor drug target, and CPR by a multiresolution computational approach. We find that upon binding to CPR, the CYP 1A1 catalytic domain becomes less embedded in the membrane and reorients, indicating that CPR may affect ligand passage to the CYP active site. Despite the constraints imposed by membrane binding, we identify several arrangements of CPR around CYP 1A1 that are compatible with ET. In the complexes, the interactions of the CPR FMN domain with the proximal side of CYP 1A1 are supplemented by more transient interactions of the CPR NADP domain with the distal side of CYP 1A1. Computed ET rates and pathways agree well with available experimental data and suggest why the CYP-CPR ET rates are low compared to those of soluble bacterial CYPs.

The Effect of Force-Field Parameters on Cytochrome P450-Membrane Interactions: Structure and Dynamics.

The simulation of membrane proteins requires compatible protein and lipid force fields that reproduce the properties of both the protein and the lipid bilayer. Cytochrome P450 enzymes are bitopic membrane proteins with a transmembrane helical anchor and a large cytosolic globular domain that dips into the membrane. As such, they are representative and challenging examples of membrane proteins for simulations, displaying features of both peripheral and integral membrane proteins. We performed molecular dynamics simulations of three cytochrome P450 isoforms (2C9, 2C19 and 1A1) in a 2-oleoyl-1-palmitoyl-sn-glycerol-3-phosphocholine bilayer using two AMBER force field combinations: GAFF-LIPID with ff99SB for the protein, and LIPID14 with ff14SB for the protein. Comparison of the structural and dynamic properties of the proteins, the lipids and the protein-membrane interactions shows differing sensitivity of the cytochrome P450 isoforms to the choice of force field, with generally better agreement with experiment for the LIPID14 + ff14SB combination.

Differing Membrane Interactions of Two Highly Similar Drug-Metabolizing Cytochrome P450 Isoforms: CYP 2C9 and CYP 2C19.

The human cytochrome P450 (CYP) 2C9 and 2C19 enzymes are two highly similar isoforms with key roles in drug metabolism. They are anchored to the endoplasmic reticulum membrane by their N-terminal transmembrane helix and interactions of their cytoplasmic globular domain with the membrane. However, their crystal structures were determined after N-terminal truncation and mutating residues in the globular domain that contact the membrane. Therefore, the CYP-membrane interactions are not structurally well-characterized and their dynamics and the influence of membrane interactions on CYP function are not well understood. We describe herein the modeling and simulation of CYP 2C9 and CYP 2C19 in a phospholipid bilayer. The simulations revealed that, despite high sequence conservation, the small sequence and structural differences between the two isoforms altered the interactions and orientations of the CYPs in the membrane bilayer. We identified residues (including K72, P73, and I99 in CYP 2C9 and E72, R73, and H99 in CYP 2C19) at the protein-membrane interface that contribute not only to the differing orientations adopted by the two isoforms in the membrane, but also to their differing substrate specificities by affecting the substrate access tunnels. Our findings provide a mechanistic interpretation of experimentally observed effects of mutagenesis on substrate selectivity.

Influence of Transmembrane Helix Mutations on Cytochrome P450-Membrane Interactions and Function.

Human cytochrome P450 (CYP) enzymes play an important role in the metabolism of drugs, steroids, fatty acids, and xenobiotics. Microsomal CYPs are anchored in the endoplasmic reticulum membrane by an N-terminal transmembrane (TM) helix that is connected to the globular catalytic domain by a flexible linker sequence. However, the structural and functional importance of the TM-helix is unclear because it has been shown that CYPs can still associate with the membrane and have enzymatic activity in reconstituted systems after truncation or modification of the N-terminal sequence. Here, we investigated the effect of mutations in the N-terminal TM-helix residues of two human steroidogenic enzymes, CYP 17A1 and CYP 19A1, that are major drug targets for cancer therapy. These mutations were originally introduced to increase the expression of the proteins in Escherichia coli. To investigate the effect of the mutations on protein-membrane interactions and function, we carried out coarse-grained and all-atom molecular dynamics simulations of the CYPs in a phospholipid bilayer. We confirmed the orientations of the globular domain in the membrane observed in the simulations by linear dichroism measurements in a Nanodisc. Whereas the behavior of CYP 19A1 was rather insensitive to truncation of the TM-helix, mutations in the TM-helix of CYP 17A1, especially W2A and E3L, led to a gradual drifting of the TM-helix out of the hydrophobic core of the membrane. This instability of the TM-helix could affect interactions with the allosteric redox partner, cytochrome b5, required for CYP 17A1's lyase activity. Furthermore, the simulations showed that the mutant TM-helix influenced the membrane interactions of the CYP 17A1 globular domain. In some simulations, the mutated TM-helix obstructed the substrate access tunnel from the membrane to the CYP active site, indicating a possible effect on enzyme function.

Inter-subunit interactions drive divergent dynamics in mammalian and Plasmodium actin filaments.

Cell motility is essential for protozoan and metazoan organisms and typically relies on the dynamic turnover of actin filaments. In metazoans, monomeric actin polymerises into usually long and stable filaments, while some protozoans form only short and highly dynamic actin filaments. These different dynamics are partly due to the different sets of actin regulatory proteins and partly due to the sequence of actin itself. Here we probe the interactions of actin subunits within divergent actin filaments using a comparative dynamic molecular model and explore their functions using Plasmodium, the protozoan causing malaria, and mouse melanoma derived B16-F1 cells as model systems. Parasite actin tagged to a fluorescent protein (FP) did not incorporate into mammalian actin filaments, and rabbit actin-FP did not incorporate into parasite actin filaments. However, exchanging the most divergent region of actin subdomain 3 allowed such reciprocal incorporation. The exchange of a single amino acid residue in subdomain 2 (N41H) of Plasmodium actin markedly improved incorporation into mammalian filaments. In the parasite, modification of most subunit-subunit interaction sites was lethal, whereas changes in actin subdomains 1 and 4 reduced efficient parasite motility and hence mosquito organ penetration. The strong penetration defects could be rescued by overexpression of the actin filament regulator coronin. Through these comparative approaches we identified an essential and common contributor, subdomain 3, which drives the differential dynamic behaviour of two highly divergent eukaryotic actins in motile cells.

An evaluation of the CYP2D6 and CYP3A4 inhibition potential of metoprolol metabolites and their contribution to drug-drug and drug-herb interaction by LC-ESI/MS/MS.

The aim of the present study was to evaluate the contribution of metabolites to drug-drug interaction and drug-herb interaction using the inhibition of CYP2D6 and CYP3A4 by metoprolol (MET) and its metabolites. The peak concentrations of unbound plasma concentration of MET, α-hydroxy metoprolol (HM), O-desmethyl metoprolol (ODM) and N-desisopropyl metoprolol (DIM) were 90.37 ± 2.69, 33.32 ± 1.92, 16.93 ± 1.70 and 7.96 ± 0.94 ng/mL, respectively. The metabolites identified, HM and ODM, had a ratio of metabolic area under the concentration-time curve (AUC) to parent AUC of ≥0.25 when either total or unbound concentration of metabolite was considered. In vitro CYP2D6 and CYP3A4 inhibition by MET, HM and ODM study revealed that MET, HM and ODM were not inhibitors of CYP3A4-catalyzed midazolam metabolism and CYP2D6-catalyzed dextromethorphan metabolism. However, DIM only met the criteria of >10% of the total drug related material and <25% of the parent using unbound concentrations. If CYP inhibition testing is solely based on metabolite exposure, DIM metabolite would probably not be considered. However, the present study has demonstrated that DIM contributes significantly to in vitro drug-drug interaction. Copyright © 2016 John Wiley & Sons, Ltd.

Identification of leads for antiproliferative activity on MDA-MB-435 human breast cancer cells through pharmacophore and CYP1A1-mediated metabolism.

CYP1A1 is a potential target for anticancer drug development due to its overexpression in certain cancer cells and role in cancer progression. To identify new leads for CYP1A1 mediated anticancer action, we attempted ligand based pharmacophore mapping, virtual screening of databases, molecular docking, MetaSite based filtering, and molecular dynamics simulations. Initial computational and in vitro screening identified 11 compounds from which we identified two lead compounds, ZINC33468944 and ZINC32101539, showed potential antitumor activity on MDA-MB-435 cell lines (GI50 < 0.1 µM) and CYP1A1 inhibition of 0.13 and 0.3 µM, respectively. Furthermore, the lead compounds were evaluated for CYP1A1 mediated metabolism, showing N-hydroxylated metabolites, which have potential of DNA adduct formation and cause cancerous cell death. Analysis of molecular dynamics simulations provided important guidelines for the further modification of the lead compounds. Hence, we claim the lead molecules for further development in anticancer drug discovery.

Synthesis and biology of ring-modified l-Histidine containing thyrotropin-releasing hormone (TRH) analogues.

Thyrotropin-releasing hormone (TRH) analogues bearing halogen groups (Cl, Br and I) at the C-2 and/or C-5 position, and the alkyl group (CH3, C2H5, C3H7, CH2C6H5) at the N-1 position of the imidazole ring of the central histidine residue were synthesized and evaluated for the receptor binding, calcium mobilization (FLIPR), and IP-1 assay at the HEK mTRHR1 and HEK mTRHR2 expressing cell lines. The most promising analogue 7k showed 925-fold selectivity for HEK mTRH-R2 receptor subtype in the IP-1 assay, 272-fold selectivity for HEK mTRH-R2 receptor subtype in the FLIPR assay, and 21-fold receptor binding specificity at HEK TRH-R2 receptor subtype. The peptide 7k was evaluated in vitro in a brain membrane competitive binding assay, and for stability analysis in the presence of TRH-DE, in vivo. The analogue 7k showed decrease in the sleeping time by more than 76% in a pentobarbital-induced sleeping assay, and showed comparatively less elevation in the TSH level in the blood, in vivo. The computational homology modeling of TRH-R1 and TRH-R2 and docking study with the most potent peptide 7k provide impetus to design CNS specific TRH analogues.

Characterization of differences in substrate specificity among CYP1A1, CYP1A2 and CYP1B1: an integrated approach employing molecular docking and molecular dynamics simulations.

Recent trends in new drug discovery of anticancer drugs have made oncologists more aware of the fact that the new drug discovery must target the developing mechanism of tumorigenesis to improve the therapeutic efficacy of antineoplastic drugs. The drugs designed are expected to have high affinity towards the novel targets selectively. Current research highlights overexpression of CYP450s, particularly cytochrome P450 1A1 (CYP1A1), in tumour cells, representing a novel target for anticancer therapy. However, the CYP1 family is identified as posing significant problems in selectivity of anticancer molecules towards CYP1A1. Three members have been identified in the human CYP1 family: CYP1A1, CYP1A2 and CYP1B1. Although sequences of the three isoform have high sequence identity, they have distinct substrate specificities. The understanding of macromolecular features that govern substrate specificity is required to understand the interplay between the protein function and dynamics, design novel antitumour compounds that could be specifically metabolized by only CYP1A1 to mediate their antitumour activity and elucidate the reasons for differences in substrate specificity profile among the three proteins. In the present study, we employed a combination of computational methodologies: molecular docking and molecular dynamics simulations. We utilized eight substrates for elucidating the difference in substrate specificity of the three isoforms. Lastly, we conclude that the substrate specificity of a particular substrate depends upon the type of the active site residues, the dynamic motions in the protein structure upon ligand binding and the physico-chemical characteristics of a particular ligand. Copyright © 2016 John Wiley & Sons, Ltd.

Ligand tunnels in T. brucei and human CYP51: Insights for parasite-specific drug design.

BACKGROUND: Cytochrome P450 sterol 14α-demethylase (CYP51) is an essential enzyme for sterol biosynthesis and a target for anti-parasitic drug design. However, the design of parasite-specific drugs that inhibit parasitic CYP51 without severe side effects remains challenging. The active site of CYP51 is situated in the interior of the protein. Here, we characterize the potential ligand egress routes and mechanisms in Trypanosoma brucei and human CYP51 enzymes. METHODS: We performed Random Acceleration Molecular Dynamics simulations of the egress of four different ligands from the active site of models of soluble and membrane-bound T. brucei CYP51 and of soluble human CYP51. RESULTS: In the simulations, tunnel 2f, which leads to the membrane, was found to be the predominant ligand egress tunnel for all the ligands studied. Tunnels S, 1 and W, which lead to the cytosol, were also used in T. brucei CYP51, whereas tunnel 1 was the only other tunnel used significantly in human CYP51. The common tunnels found previously in other CYPs were barely used. The ligand egress times were shorter for human than T. brucei CYP51, suggesting lower barriers to ligand passage. Two gating residues, F105 and M460, in T. brucei CYP51 that modulate the opening of tunnels 2f and S were identified. CONCLUSIONS: Although the main egress tunnel was the same, differences in the tunnel-lining residues, ligand passage and tunnel usage were found between T. brucei and human CYP51s. GENERAL SIGNIFICANCE: The results provide a basis for the design of selective anti-parasitic agents targeting the ligand tunnels.

Synthesis of CNS active thyrotropin-releasing hormone (TRH)-like peptides: Biological evaluation and effect on cognitive impairment induced by cerebral ischemia in mice.

Thyrotropin-releasing hormone (TRH)-like peptides were synthesized by replacing critical histidine and pGlu residues in the native peptide. The peptides were evaluated in vitro for receptor binding activity assay and in the cell functional assay; the peptides exhibit selective basal signaling agonist behavior toward TRH-R2. For example, peptides 8a, 8b, 8c, 8 f, 8 h, 8 l and 12 d activated TRH-R2 with potency (EC50) of 0.53 µM, 0.048 µM, 0.05 µM, 0.006 µM, 0.31 µM, 0.034 µM and 0.004 µM, respectively. In contrast for signaling activation of TRH-R1, the same peptide required higher concentration of 19.35 µM, 3.98 µM, 2.54 µM, 0.287 µM, 11.28 µM, 0.986 µM and 0.944 µM, respectively. The results showed that peptides were 36.5, 82.9, 50.8, 47.8, 36.3, 32.6 and 235-fold selective to TRH-R2 receptor subtype. The peptides were investigated for CNS activity at 10 µmol/kg in pentobarbital-induced sleep assay study. Peptides 8c (16.5 ± 1.4 min) and 8l (16.5 ± 2.1 min) displayed excellent CNS activity. In an in vivo study, peptide 8c did not cause significant change in the rat plasma TSH levels. The peptide 8c was further investigated for neuroprotective potential, and significantly reduced infracts volume and neurological score in the focal cerebral ischemia model in mice. Peptide 8c also significantly lowered MDA levels, indicating reduction of oxidative and enhanced percentage cell survival in CA1 region, when compared to ischemic brain.

Plasma protein binding, pharmacokinetics, tissue distribution and CYP450 biotransformation studies of fidarestat by ultra high performance liquid chromatography-high resolution mass spectrometry.

Fidarestat, an aldose reductase inhibitor, has been used for the treatment of the diabetic associated complications such as retinopathy, neuropathy and nephropathy. To better understand the metabolism and pharmacokinetics of fidarestat, we have evaluated plasma protein binding, pharmacokinetics, tissue distribution of the drug and its conjugated metabolites and CYP450 biotransformation by liquid chromatography-high resolution mass spectrometry. Effective chromatographic separation of fidarestat and hydrochlorothiazide (IS) in rat plasma and tissues was achieved on Hypersil gold C-18 column in an isocratic elution mode. For detection, a high-resolution Orbitrap mass spectrometer with heated electrospray ionization inlet in the negative ion mode was used. High-resolution extracted ion chromatograms for each analyte were obtained by processing the full-scan MS mode with 5 ppm mass tolerance. The impact of plasma protein binding with the drug and conjugated metabolites of the drug on pharmacokinetics has been determined. The study indicated that 9.5% of free form of fidarestat may be pharmacologically active and the Cmax for free fidarestat was found to be 80.30 ± 6.78 ng/mL. The AUC0-t and AUC0-∞ were found to be 185.46 ± 32 and 195.92 ± 15.06 ng h/mL, respectively. Among tissues, the maximum observed distribution was found to be in kidney followed by liver and heart. Docking experiments and in vitro CYP450 reaction phenotyping revealed that two CYP1A2 and CYP2D6 are involved in the phase I metabolism of fidarestat. Oxidative deamination and N/O glucuronidation are the major phase I and phase II metabolites, respectively. In vitro CYP450 inhibition assay of fidarestat for drug-drug interaction showed weak inhibition and may not alter pharmacokinetics, distribution or clearance of other co-administered drug.

On the application of the MARTINI coarse-grained model to immersion of a protein in a phospholipid bilayer.

An important step in the simulation of a membrane protein in a phospholipid bilayer is the correct immersion of the protein in the bilayer. Crystal structures are determined without the bilayer. Particularly for proteins with monotopic domains, it can be unclear how deeply and in which orientation the protein is being inserted in the membrane. We have previously developed a procedure combining coarse-grain (CG) with all-atom (AA) molecular dynamics (MD) simulations to insert and simulate a cytochrome P450 (CYP) possessing an N-terminal transmembrane helix connected by a flexible linker region to a globular domain that dips into the membrane. The CG simulations provide a computationally efficient means to explore different orientations and conformations of the CYP in the membrane. Converged configurations obtained in the CG simulations are then refined in AA simulations. Here, we tested different variants of the MARTINI CG model, differing in the water model, the treatment of long-range non-bonded interactions, and the implementation (GROMACS 4.5.5 vs 5.0.4), for this purpose. We examined the behavior of the models for simulating a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer in water and for the immersion of CYP3A4 in a POPC bilayer, and compared the CG-MD results with the previously reported experimental and simulation results. We also tested the methodology on a set of four other CYPs. Finally, we propose an optimized protocol for modeling such protein-membrane systems that provides the most plausible configurations and is computationally efficient; this incorporates the standard non-polar water model and the GROMACS 5.0.4 implementation with a reaction field treatment of long-range interactions.

Comparative proteomics among cytochrome p450 family 1 for differential substrate specificity.

Apart from playing key roles in drug metabolism and adverse drug-drug interactions, CYPs are potential drug targets to treat a variety of diseases. The intervention of over expression of P450 1A1 (CYP1A1) in tumor cells is identified as a novel strategy for anticancer therapy. We investigated three isoforms of CYP1 family (CYP1A1, CYP1A2, and CYP1B1) for their substrate specificity. The understanding of macromolecular features that govern substrate specificity is required to understand the interplay between the protein function and dynamics. This can help in design of new antitumor molecule specifically metabolized by CYP1A1 to mediate their antitumor activity. In the present study, we carried out the comparative protein structure analysis of the three isoforms. Sequence alignment, root mean square deviation (RMSD) analysis, B-factor analysis was performed to give a better understanding of the macromolecular features involved in substrate specificity and to understand the interplay between protein dynamics and functions which will have important implications on rational design of anticancer drugs. We identified the differences in amino acid residues among the three isoforms of CYP1 family, which may account for differential substrate specificity. Six putative substrate recognition sequences are characterized along with the regions they form in the protein structure. Further the RMSD and B-factor analysis provides the information about the identified residues having the maximum RMSD and B-factor deviations.

Pregnane X Receptor and P-glycoprotein: a connexion for Alzheimer's disease management.

The translational failure between preclinical animal models and clinical outcome has alarmed us to search for a new strategy in the treatment of Alzheimer's disease (AD). Interlink between Pregnane X Receptor (PXR) and P-glycoprotein (Pgp) at the blood brain barrier (BBB) has raised hope toward a new disease modifying therapy in AD. Pgp is a major efflux transporter for beta amyloid (Aß) at human BBB. A literature survey reveals diminished expression of Pgp transporter at the BBB in AD patients. Pregnane X Receptor is a major transcriptional regulator of Pgp. Restoration of Pgp at the BBB enhances the elimination of the Aß from brain alongside and inhibits the apical to basolateral movement of Aß from the circulatory blood. This review concentrates on in vitro, in vivo, and in silico advancements on the study of the PXR in context to Pgp and discusses the substrate and inhibitor specificity between PXR and Pgp.

Intestinal transport of TRH analogs through PepT1: the role of in silico and in vitro modeling.

The present study involves molecular docking, molecular dynamics (MD) simulation studies, and Caco-2 cell monolayer permeability assay to investigate the effect of structural modifications on PepT1-mediated transport of thyrotropin releasing hormone (TRH) analogs. Molecular docking of four TRH analogs was performed using a homology model of human PepT1 followed by subsequent MD simulation studies. Caco-2 cell monolayer permeability studies of four TRH analogs were performed at apical to basolateral and basolateral to apical directions. Inhibition experiments were carried out using Gly-Sar, a typical PepT1 substrate, to confirm the PepT1-mediated transport mechanism of TRH analogs. Papp of the four analogs follows the order: NP-1894 < NP-2378 < NP-1896 < NP-1895. Higher absorptive transport was observed in the case of TRH analogs, indicating the possibility of a carrier-mediated transport mechanism. Further, the significant inhibition of the uptake of Gly-Sar by TRH analogs confirmed the PepT1-mediated transport mechanism. Glide docking scores of all the four analogues were in good agreement with their transport rates, suggesting the role of substrate binding affinity in the PepT1-mediated transport of TRH analogs. MD simulation studies revealed that the polar interactions with amino acid residues present in the active site are primarily responsible for substrate binding, and a downward trend was observed with the increase in bulkiness at the N-histidyl moiety of TRH analogs.

Predicting drug metabolism by CYP1A1, CYP1A2, and CYP1B1: insights from MetaSite, molecular docking and quantum chemical calculations.

Recently, CYP1 enzymes are documented for selective metabolism of anticancer leads in cancer prevention and/or progression. Elucidation of specificity of substrates/inhibitors of CYP1 isoforms plays a vital role in design of more selective and potent anticancer leads. However, an area of concern is the broad range of substrate specificities and planar nature of substrates with limited dataset which makes it difficult to predict their site of metabolism (SOM) accurately. In the present study, various models for prediction of site of metabolism in case of CYP1A1, CYP1A2, and CYP1B1 substrates were developed using MetaSite, molecular docking, and quantum chemical descriptors. The predictive accuracy of MetaSite, molecular docking, and quantum chemical descriptors in identifying experimental site of metabolism was analyzed at three levels; top rank, top three ranks, and top five ranks. Two quantum chemical descriptors, chemical hardness and local nucleophilicity are proposed for the prediction of CYP-mediated SOM for the first time. The predictive accuracy shown by chemical hardness at top three ranks was 83.3, 85.7, and 84.6 % for CYP1A1, CYP1A2 and CYP1B1, respectively, whereas local nucleophilicity gave poor predictions of 50, 42.8, and 46.2 %, respectively. The predictability of chemical hardness descriptor outperformed at all three levels of ranks for CYP1A1, CYP1A2, and CYP1B1. Hence, we propose chemical hardness as an useful quantum chemical descriptor for prediction of metabolically vulnerable prints in CYP1A1, CYP1A2, and CYP1B1 mediated metabolism and support the optimization efforts in drug discovery and development programs.